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The Polchinski version of the exact renomiaUzation group equation is discussed and 
its appHcations in scalar and ferniionic theories are reviewed. Relation between 
this approach and the standard renormalization group is studied, in particular the 
relation between the derivative expansion and the perturbation theory expansion 
is worked out in some detail. 

C^^ ' 1 Introduction 

o\\ 

O^ ■ As it is known (see an interesting review by Shirkovl) thejiasic idea of the renormalization group (RG) was 

\ formulated in an article by Stueckelberg and PetermannEI. The existence of such group of transformations 

i-O ■ was related to the arbitrariness in the procedure of subtraction of ultraviolet divergencies in quantum 

r-y I electrodynamics. Functional equations for the propagators in the ultraviolet limit, coxresponding to 

^"^ ■ the regularized Dyson transformations, were basically derived by Gell-Mann and LowH. Bogoliubov 

^v| \ and ShirkovQ unveiled the group nature of these functional equations, established their relation to the 

■ renormalization group by Stueckelberg and Petermann and derived the functional group equations for 

*^ I propagators and vertices of quantum electrodynamics for a general case, i.e. massive case. For example, 

^ ■ in a renormalizable quantum field theory model with one coupling constant g and one mass m (like gc/)'* 

04 \ model in four dimensions) the functional equation for the invariant charge g{x, y; g), with x = p^ / ^^ and 

2~ ■ y — m? I ^^ being the ratios of the momentum and mass squared to the square of the normalization scale, 

~^ ' has the following form: 

vN ' ^ X V \ 

O; 9{x,y]g)=9y- , - , g{t,y]g)j . (1) 

^^ , This equation establishes the connection between the exact symmetry of Stueckelberg and Petermann 

^^ ' and an approximate one of Gell-Mann and Low since in the massless limit to = it reduces formally to 

-4— > , the result by Gell-Mann and LowH. Eq. ([l|) reflects the exact symmetry of a solution which later became 

O ' known as a functional self-similarity symmetry Q (as it is relevant to the notion of self-similarity known 

O . in mathematical physics). From this equation the standard differential renormalization group equations 

i-Q ' in both Gell-Mann-Low and Callan-Symanzik formErll can be derived. We will call this approach just 

^ , " renormalization group" (RG) in order to distinguish it from the approach originated from the works by 

Wilson and later by Polchinski, which is traditionally called "exact renormalization group" (ERG). 

The symmetry, underlying the RG, is the symmetry of a characteristic function i^(x;xo,-fo) of a 
C^ ' physical problem ('Isolution of the problem") with respect to the change of the boundary condition 

F{xq\xq,Fq) = -Fo El. If we now consider _F(a;; a;i,i^i), which is the same solution but with boundary 
condition settled at another point xi, we arrive to the functional relation 

F(x;a;i,i^o) = F{x]Xi,F{xi\Xq,Fq)). 

Assuming that the function possesses the homogeneity property _F(a;; xq, -fo) = / ( ~i ^o ) with /(I, _Fo) = 
Fq we obtain for the function / an equation of the type dl^) . _ _ 

Another approach was developed in articles by Kadanoff and Wilson I3'U. It employes the idea of 
Wilson effective action, which is the action obtained by integrating out degrees of freedom with momenta 
p2 ^ ^2 jj^ ^j^g defining functional integral. So, when the scale A is reduced, the generating functional 



X 



ZA(J) = y n|p|<Ai?,^(p)exp|-5(</>;A)-y dp J<p\ (2) 

includes integration over less and less modes, but the effective action S'(0; A) changes in such a way that 
the generating functional remains unchanged. In this sense, it describes the same physics for any A, that 
is the ^-matrix elements and Green functions remain unchanged. These ideas were first developed for 
lattice statistical systems and were fruitfully used in condensed matter physics. Some re-formulations of 
this approach, which made possible its applications in quantum field theory, were made first by Wegner 



and HoughtonE3 and WeinbergEJ and later by PolchinskiE3. It is basically this approach which was called 
the exact renormalization group (ERG). In this article_sve discuss in detail the Polchinski version of the 
ERG (see a comprehensive review by Ball and ThorneE^ on this subject), though results obtained within 
other schemes will also be presented. 

The action can be characterized by a set g — {32, <?4, • ■ •} of coupling constants of all possible oper- 
ators consistent with the symmetries of the system. The change of the scale A -^ A/t combined with 
the corresponding change of the set g — > Rtg, where Rt is some operator, can be viewed as a group 
transformation. Under such transformation the field changes as (/)(p) —>■ at4i{tp). Then, for example, for 
the 2-point function, defined from (g) in a standard way, we have a relation of the type: 

G2(|,ff(A))=a?G2(^|,i?*<?(A)). (3) 

One can see similarity of this relation with the functional equation (nl). Having in mind this similarity, 
we can view the symmetry, underlying the Wilson renormalization group, as the invariancc of the Wilson 
effective action with respect to the boundary condition 

5(0,A)|a=Ao=^o(</'), (4) 

where So{4>) is some "fundamental" action defined at the "fundamental" scale Aq. Because of this 
similarity one can expect that the RG and ERG approaches have much in common. This aspect was 
emphasized by ShirkovEJ. 

The first aim of this contribution is to give a short review of recent applications of the ERG in quantmn 
field theory. Here we will limit ourselves to the case of scalar and fermionic theories only. Application 
of the ERG in gauge theories at the moment encounters difficulties (which are not unavoidable, in our 
opinion) due to the absence of a gauge invariant formulationlij. We will concentrate on results which are 
non-perturbative, though obtained within some approximation. This approximation will be the derivative 
expansion. 

The second aim is to work out some details on the relation between the RG and the ERG. To much 
extent this was studied in the original paper by Polchinskili3 (for a detailed discussion see also the article 
by Ball and ThorneE3) . Explicit calculation of the 1-loop /^function and anomalous dimensions in scalar 
theory with the ERG method was done by Hughes and LiuE3. In the present contribution we will examine 
the /3-function calculated by the derivative expansion technique within the ERG and analyze its relation 
to the perturbative RG calculation. 

The article is organized as follows. In Sect. 2 we review the Polchinski version of the ERG equation 
in the case of scalar theory and describe the derivative expansion, which is a method used for non- 
perturbative calculations in this approach. In Sect. 3 we explain how the standard perturbative results 
can be obtained in the framework of the ERG approach and compare them with the results calculated 
within the derivative expansion. Wc will discuss the difference and the relation between these two 
calculations. In Sect. 4 we review very briefly the main non-perturbative results obtained in scalar field 
theory within the ERG approach. Here calculations for various types of the ERG equations, not only for 
the Polchinski one, are summarized. In the last section the ERG equation for fermions and some first 
results obtained in the 2-dimensional Gross-Neveu type model are presented. 

2 The ERG equation 

Consider a d-dimensional scalar theory of the field $ which is invariant ujider the Z2-symmetry transfor- 
mation $ ^ — $. Let us introduce the regularized generating functional t5ll3 



^-W{J;A) 



I D^e-^^^'^'-'\ (5) 



where the action with the source term is given by 

Jd 



d^pJ{p)Q],\pMp) + fA. (6) 



The regularized propagator 



Pa{p) 



K 



m 



pz 



(7) 



is defined by introducing the regulating function K{z) which is supposed to be decreasing fast enough 
when z — > cx) and is normalized by K{0) — 1. The functions Qk{p) and /a are necessary for the 
consistency of the formulation and are to be determined. Note that the last term in Eq. (pi) does not 
depend on the field. 

In the functional integral above only those momentum modes which have \p\ lower or about A are 
important, contributions of higher modes are suppressed by the the regulating function. Thus A plays the 
role of a (smooth) upper cutoff. In this setting Wilson's idea is realized as follows: the effective action S'a 
is such that while the scale A varies the S'-matrix elements or even the off-shell Green functions remain 
unchanged. This implies that the generating functional is a constant function of A: 



A 



dW(J,K) 
dK 



0. 



(8) 



The change of the propagator (the range of modes suppressed in the functional integral) with A in the 
kinetic term is compensated by the change of the action of interaction and other terms in Eq. (0) so that 
the whole functional integral, defining the generating functional, describes the same physics. Next step 
is to use the functional integral identity 



D^- 



5^{p) V2(5$(-p) 



Pl'{p)'^{p) + Ql\p)J{p) e-^-(*'^) = 0. 



Combining this identity with the condition (pi) one can derive the ERG equation for the effective action. 
Here we omit details of the derivation and present the result only. Before doing this let us introduce some 
notations. It is convenient to define the parameter 



t 



In 



A 
A^' 



(9) 



where Aq is some fixed scale. The scale of dimensionful objects is carried by A. So we can define the 
" dimensionless momentum" q and the dimensionlcss field variable 4>{q]t) as follows: 



P_ 

A' 



;t) = Ai+'^/2$(p;A,Ao). 



The dependence of the field on t is characterized by the anomalous dimension 77: 



|0(.) = Im. 



(10) 



(11) 



Then the part of the effective action 
/a, takes the form 



which does not include the source term and the constant term 



S\^ 



d'^q q 

d'^(9) 



(2^) 



The ERG equation for it is 
'dt 



d'^q{2TTfK'{q^ 
SS 



K{q^) 



(f){-q) + Sint[(f>;t\. 



SS SS 



d^q^biq) 



H{q) 



S(t>{q)S^{-q) S^{q) S(P{-q) 

K'jq') d 

K{q^) ^ 2 



Sd 



(12) 



(13) 



2q 



d'^q(t>{q)q,. 



d' SS 

9?^ ^<Piq) ' 



The prime in the last term means that the derivative does not act on the standard (5-functions of the total 
energy-momentum conservation which appear in the action in the momentum representation. Only the 
first line is a non-trivial part, the rest of the terms just refiects the canonical dimensions of the objects 
of the action and the anomalous dimension of the field. For the condition ^ to fulfill the function /a 
has to satisfy 



fA = Jd^'p 



2f 2^ PA{p) 

Plip) 



Q-'ip') 



J{p)J{-p) - I d'pj^SiO) 



and Qa{p) — Pa[p)Q{p^)i where Q{p^) obeys the equation Q{p^) = [ri/2)Q{p'^). Here the dot means 
differentiation with respect to t. 

The equation dl3) is supposed to be suppUed with the initial condition set at some scale Aq or at 

Sint[4'\i]\t=0 = Sint[4>] = / d''-xLinti(t>), 



where L^nt is essentially the bare Lagrangian. Then the ERG equation defines the running Lagrangian 
Lint{4>T't) and, correspondingly, the running action Sint[4>,t] = ^ cfixLint{4>,t), i-e. a trajectory in the 
space of all possible Lagrangians parametrized by t. Note that the equation (|l^) is exact and non- 
perturbative. One can observe certain similarity of it with a functional generalization of tJie RG equation 
to the case of Lagrangians of arbitrary type (including non-renormalizable Lagrangians) Il3. 

The limit t ^ oo describes either the situation with Aq fixed and A — > 0, i.e. the limit of low 
characteristic energy, or the situation Aq -^ oo, i.e. the continuum limit of the model. In this way the 
ERG equation allows for non-perturbative studies of the continuum limit in quantum field theory. 

There are a few important issues that can be addressed by studying and solving this equation. First of 
all we can look for fixed point solutions L*^^{(j)) which some of the trajectories can approach as i — > +cxd 
(here we mean fixed points for finite values of the coupling constants in the Lagrangian). The fixed 
points satisfy Eq. (03) with the zero l.h.s., S'*^^ — 0, and this equation defines the value of the anomalous 
dimension 77 = if for which such solution exists. The Gaussian fixed point 5*„j = gives an example of 
the trivial solution. 

Having found a fixed point solution wc can study the theory in its vicinity. For this we represent the 
Lagrangian as an expansion in operators 



i^„t(</>,t)=i*„t('/>) + > One^"^ (14) 






The parameters A„ are called critical exponents and they are physical observables. The operators On (</)), 
which correspond to A„ > 0, are called relevant operators and are important for the physics of the system 
in the vicinity of the fixed point when t —^ 00. The ERG equation ( [l3| ) in principle allows to calculate 
the critical exponents and found the corresponding operators. 

Finally, one can try to solve the ERG equation for arbitrary i, i.e. find the complete renormalization 
group trajectory. There are also other interesting problems which can be considered in a non - perturbative 
way in the framework of this approach. These include bound states, Zamolodchikov c-function, etc. 

We would like to mention that in quantum field theory ERG equaiipns of other types are also 
considered. Historically the first one was the Wegner- Houghton equation 113. It was formulated for the 
sharp cutoff and was used in a number of articles for calculation of fixed points, critical exponents, 
flows, etc. in the scalar theory (see the article by Hasenfratz and Hasenfratz for one of the first detailed 
studies tZi). An approach based on an ERG equation with a sharp cutoff for the cfFective_action T^ff was 
developed by MorrisEfl. Another version is the equation for the average effective actionE£l. Some of these 
results will be discussed in Sect. 4. 

3 Approximations and relation between ERG and RG 

There are some special cases when the ERG equation ( p^ ) simplifies essentially and one can find its 
iLujtions. This happens, for example, in the theory with iV-component scalar field in the limit N —>■ 00 



H23. In general it is not clear how the ERG equation (13) can be solved exactly. So we need to use 
an approximation to obtain solutions and to analyze them. For this it is useful to have an idea of a 
Lagrangian which captures essential features of the problem. In general, even if we start with a simple 
initial Lagrangian, like, for example, Lint{4>) = g4>'^/4\ at i = 0, the running Lagrangian can include 
all possible operators constructed out of the field 4> and its derivatives which are consistent with the 
symmetry of the problem. In the momentum representation the action can be written as an infinite series 

Srnt[(b,t] = -J.^—^A2iq,~q,t)q^{q)<j>i~q) + -J ^^^^ Si}_^q,) 

X A4{qi,q2,q3,q4,t)c/){qi)c/){q2)(f>{q3)(l)(q4) + ■■■ (15) 

4 



Then the ERG equation gives rise to an infinite system of coupled equations 

A2{q,-q,t) = [2 + ri)A2{q,-q,t)-2q^^—A2[q,-q,t) + 2K'{q^)A2{q,-q,t) 



dqn 



d'li r.> 



{2n) 



-^K'{q-^)Ai{q, -q, qi, -qi,t), 



d 



A^{qi,...,q^,t) = (4-d + 2r7)A4-^qj>— — A4 

4 

+ 2^ [K'{q])A2{qj, -q,,t)\ ^4(^1, . . . 94, t) 



(2^) 



-^K'{q )Ae{qi,...,q4,q,-q,t), 



(16) 



Let us consider now the case d = 4 as an example. The initial conditions, corresponding to one of the 
jil renormalization prescriptions in terms of the coefficient functions in (15), are settled at two scales 



at A = Aq: 



A2{q,-q,Ao) = pi(Ao) + gV2(Ao); 
-44(91, ■•■, 94, Ao) = 5(Ao) = .gs; 
A2j{qi, • . • , g2j, Ao) = for j > 3, 



(17) 



that gives the standard form for the bare Lagrangian where, of course, gs is the bare coupling constant; 
and at some physically relevant scale (renormalization point) A = /Xi?: 



^4(0,0,0,0,/ifl) 
Pi (Mi?) = 0, 



5(Mfl); 

P2(Mi?) = 0. 



(18) 



In the formulas above for illustrative purposes we indicated the dependence on the scale as the dependence 
on A and not on t as before. 

One of the known ways to solve the ERG equation is to use the perturbative expansion. We assume 
that all coefhcient functions can be represented as power series in gR = g{ij,fi). After a bit lengthy but 
straightforward calculation one can get solutions for A2j. In particular. 



5(A;Mfl,5fl) 



A4{Q,0,0,0,A;fiR,gR) 


' d (p^\ 

dA' VA'2y 

(t)-'<' 

)] + 0{g%). 


) 

fp' 


K 

\\ 
?// 


rp' 


'^ ^'^1 (2.)4p2p2|/A y 

„ /■ d^p 1 1 ri / . 

^--^^^y(2.)v^p42r 


VA'2 



-K 



Is? 



(19) 



(20) 



In this formula we indicated the dependence on the initial condition (i.e. the renormalizatioiip^oint) 
explicitly. Taking the limit Ao ^ 00 we get the standard /3-function in the 1-loop approximationll3: 



PidB) = ^;^5(A; /ifl,, gn) 



IGtt' 



:9l + 0{gl). 



(21) 



Can we get non-perturbative solutions of the ERG equation? At the moment the available systematic 
techniques for this are the derivative expansion and its modifications. The idea is to represent the action 
(here in the coordinate representation) as 



Sint [(!>, t] 



d''x[v{4>{x),t) + {d^(^fu{4>{x),t) 



(22) 



where V, U, Hi, H2, ■ ■ ■ depend on the field but not on its derivatives. In the momentum representation 
Eq. (p3) corresponds to the expansion in powers of momenta, so one can hope that using a few first terms 
of this expansion is justifiedjf we consider effects at low momenta (there may be some complications in 
the case of the sharp cutoff El). The function V{(j),t) is the effective (local) potential of the theory. To 
solve the ERG equation ( |l3| ) approximately we truncate the derivative expansion ( P2[ ) and substitute the 
truncated action into the ERG equation thus getting equations for 
if we consider just the leading and next-to-leading terms we haveC 



re coefficient functions. For example, 



2. 



V 



V = -aV" - 2(3U" + j{V')^ +d-V{l-j--\+(l)V' 

U = -aU" + 5{V"f + A-iUV" + 27C/V - ,;t/ + f 1 - J - ^ ) 0C/' - ^ 



where the prime means differentiation with respect to the field 
the parameters characterizing the regulating function: 



(23) 
(24) 
The equation depends explicitly on 



^2^ 



K'{q^), P 



(27r)° 



K'{q^), 7 = ^'(0), 5^K"{Q). 



(25) 



This is similar to the dependence on the renormalization scheme in the RG. However, there is a problem 
of the breaking of the reparametrizatioxLmyariance by the derivative expansion that gives rise additional 
dependence on the regulating functionoO. 

One way is to solve the ajiptem (p3|), (24) numerically. Similar systems for the ERG equations of 
other types were also studied E3. These results will be discussed in the next section. 

Another way is to do further approximation, namely, to expand the functions V , U, etc. in powers 
of fields: 



F(0,t) 
[/(</>, t) 



ai(i)02 
b2m^- 



-a2(t)<^^ 



a^m"" 



(26) 



Then the system (p^), (p^) becomes the following set of flow equations t 



Ol = 


= l^- 


h ri)ai — L2a2 — 002/81 + ia^, 


0.2 = 


- (4- 


-d+ 2ri)a2 - SOas - IO63/S1 + 160^, 


as = 


= (6- 


-2d + 3ri)a3 + 24aia3 -t- I6al, 


62 = 


= (2- 


-d+ 277)62 + IGaifla + 160162 - 2O63, 


k -- 


= (4- 


192 144 
-2d + 377)63 -1- -—0262 -t- 24ai63 + 24aia3 -t- -^aj 





(27) 



with the relation 



= 7/ -52(1262 -8a?). 



(28) 

Here we introduced the combinations of the scheme parameters si = aj/pS and §2 = S/j'^. The last 
equation arises because the normalization of the kinetic term is fixed. Let us mention that Eqs. (E^) are 
basically /3-functions for the coupling constants of corresponding operators. 

The leading order of the derivative expansion with subsequent polynomial apprnYimatinn nf fVip 
potential V{(f>,t) in scalar theories for the Wegner-Houghton ERG equation was studied in detail EZlo'EJ. 
We summarize some of the results in the next section. Now we are going to illustrate the relation between 
the solutions within the derivative expansion in the ERG and the perturbative results in the RG (for 
d = 4). For this purpose let us solve Eqs. (p7|), ( p8| ) perturbatively, i.e. presenting the coefficients ai{t) 
and bi{t) and 77 as series in powers of g{t) = 4!a2(i), the coupling constant of the ^''-interaction. After a 
simple calculation up to the next-to-leading order 0{d^) of the derivative expansion one gets 



Pig) = -j^9{t) 



6{a-i+-(36 + 



0.9' + 0{g^). 



(29) 



We see that the expression (21) for the /3- function to the order g^ is not recovered within this approxi- 
mation. Moreover, the result (29h for the first coefficient of the /?- function depends on the scheme. Of 
course, what happens is that there are contributions with higher derivatives (or with higher powers of 



the momenta) in the expansion ( p2[ ) which give rise to further terms in the equations (p7|), (p8[). In the 
perturbative solution these terms provide further contributions to the expression (p9|), also to the g^ coef- 
ficient (the dots in that expression stand for these contributions) . To see the form of these contributions 
it is useful to re- write Eq. (poj), using the definitions (Eq), as follows 



/3(5) = 6<7V(0^V) 



K'{0) + \q'K"{0) 



0{g' 



Analyzing the relevant contributions of higher derivatives one can show that when they are taken into 
account the expression above sums up to 

The integral can be calculated for an arbitrary regulating function and one obtains 

P{9) = T^g' + Oig'), 

the standard result (pl|). In this calculation we used only that A'(oo) — and the normalization K{0) = 1. 

This example illustrates the relation between the derivative expansion of the ERG and the perturba- 
tion theory of the RG. On one hand the derivative expansion contains only a part of the contribution of 
the weak coupling perturbative expansion to a given order in g. It can be obtained by making the Taylor 
expansion of a part of terms in the integral over internal momenta of the corresponding Feynman dia- 
gram. Indeed, consider expression (n9h or (EU), which can be interpreted as a contribution of the 1-loop 
diagram with zero external momenta regularized according to Eq. (^). If we take the limit Aq ^ oo and 
differentiate ( [l9[) with respect to A, according to the definition (^), we obtain precisely the expression 
(|30|). Now, if for infinite Aq we expand the part in the square brackets in powers of p^ (the integral still 
remains finite), we get the expression which after differentiation with respect to A provides the result 
( |29| ) of the derivative expansion in the next-to- leading order. 

On the other hand the derivative expansion, when is not limited to the weak coupling expansion, 
contains non-perturbative information, that makes it and the ERG equation in general to be rather 
valuable tools. We will discuss the results of non-perturbative calculations in the next sections. 

4 Scalar theory: main results 

Here we discuss some of the main results of the studies within the ERG approach in the scalar theory of one 
field with Z2-symmetry (symmetry under the transformation $—>—$) on the d-dimensional Euclidean 
space. We collect results obtained by various authors and for different versions of the ERG equation. 
On general grounds different versions should give physically equivalent results, and concrete calculations 
confirm this (see the discussion of fixed point solutions and critical exponents below). However, no 
rigorous proof of the equivalence of different approaches has been given so far. 

First, there were numerous studies of fixed point solutions for various d. We discuss them in turn. 

1. Fixed points: d = A. Hasenfratz and Hasenftatz performed a study of the Wegner-Houghton 



version of the ERG equation (analog of Eq. (13)) numerically in the leading approximation C'(9*')of the 
derivative-expansion (local potential approximation) and showed that it has no non-trivial fixed point 
solutions EH. 

2. Fixed points: d — 2>. In the same article it was showed, also numerically, that in order 0{d^) 
the Wigner-Houghton equation has one non-trivial fixed point solution, the known Wilson-Fisher fixed 
point 113, which is in the universality class of the three-dimensional Ising model. This fixed point was 
also studied in the 0(9^) apacoximation, both for the Polchinski type ERG equation (fi3|)E3 and for the 
Wegner-Houghton equationEJo . 

As an example, let us review the calculation for the Polchinski type ERG equation. To the order 
0(9°) of the derivative expansion the fixed point equation is Eq. ( p3| ) without the term 2j3U" and the 
l.h.s. V = Q. It was solved numerically with two boundary conditions at (/> = 0: 

l-'(0) = 0, V"{i))=p, 

the first one just reflecting the symmetry of the theory. It was shown that a fixed point solution V*{(j)) 
regular for all finite values of exists only for a special value oi p = p* = —0.2286 ... In this approximation 



the fixed point solution y*(0) does not depend on the scheme parameters and the value of the anomalous 
dimension at the fixed point is rj* = 0. In the next-to-leading approximation 0{d^) the system of fixed 
point equations , i.e. Eqs. (p3|), ( |24| ) with zero left hand sides, was solved numerically with the following 
boundary conditions for ¥{(/)) and U{(j)) a.t (p — 0: 

V'{0) = 0, V"{0) = p, U{0) = 0, f/"(0) = 0. 

A regular non-trivial solution exists and is unique only for a special value oi p = p* and the anomalous 
dimension rj = ij* ^ 0, but now p* and rj* depend on the scheme parameters (p5[). 

3. Fixed points: 2 < d < 4. For this range of dimensions the Wegner-Houghton fixed point equation 
was analyzed in the leading order of the derivative expaiisiDn with subsequent approximation of the local 
potential V*{(f>) by a polynomial in powers of the fieldc3'EJ: 

Vj:M) = al^' + «2</'' + ■ • • + c*M(b^^'. (31) 

In this case the fixed point ERG equation reduces to a system of M algebraic equations for the coefficients 
a*. As before, this approximation gives the anomalous dimension 77* = 0, nevertheless it captures some 
general features of the multicritical fixed points below d = 4. In particular, for approximation (|3|) the 
system shows the existence of the first (M — 1) upper critical dimensions d^ = 2k/ {k — 1) = 4, 3, 8/3, . . . 
for fc = 2, 3, 4, . . . , M. Then, a,t d — dk the trivial Gaussian fixed point V^j = is a branching point of a 
new non-trivial fixed point and the Gaussian one below dk ■ 

A general feature of the polynomial approximation is the appearance of numerous spurious solutions. 
For example, for d = 3 from the numerical results, discussed above, we know that there is only one 
non-trivial fixed point, the Wilson-Fisher fixed point. But if we approximate the potential V*((/)) by a 
polynomial V^j{(j)) the system can have till M real valued non-trivial solutions for the coefficients a*, all 
of them but one being spurious. The problem of spurious solutions can be solved by analyzing solutions 
of subsequent approximations with various M and selecting the stable ones, which represent the true 
fixed points. 

The polynomial approximation permits an analytical study and reproduces some of the qualitative 
features of the structure of the fixed pointSj—However, arguments were presented which show that the 
polynomial approximation is not convergentEj, i.e. for a given cj) |V^*/(0) — T^*((/))| approaches small but 
non-zero value as M — > cx). 1—11—1 

4. Fixed points: d = 2. The case of d = 2 dimensions was studied by Morris cjO. At 0{d^) 
order of the derivative expansion he found first 10 multicritical fixed points of an infinite series which 
corresponds to unitary minimal models of conformal field theory and whose existence was conjectured by 
Zamolodchikov. 

5. Critical exponents. 

To calculate the critical exponents at a given fixed point we expand the potentials V , U, etc. around 
the fixed point solution, 

V{^,t) = V*{(t>)+SV{(t),t), f/(0,i) = f/*(0) + <5L/(0,t),..., (32) 

linearize the ERG equations (|23|), ( |24| ) with respect to the linear deviations SV, 5U , etc. and represent 
these deviations in the vicinity of the fixed point as 

<5F(0,f)=^z;„(</))e^"*. (33) 

n 

The critical exponents A„ can be found as eigenvalues of the linearized system of the ERG equations. 

Most of the study was done for critical exponents of the Wilson-Fisher fixed point for d = 3. There is 
one positive critical exponent Ai and, correspondingly, one relevant operator, the rest of A„ < (n > 2). 
We give a summary of the results of various calculations of j^ = 1/Ai , characterizing the critical exponent 
of the relevant operator, w — — A2 for the first irrelevant operator, and the anomalous dimension 77 
in Table 1. One can see that the results of different schemes and approximations are in a reasonable 
agreement. In cases when the characteristic under consideration is scheme dependent the intervals of 
values, which correspond to certain ranges of values of the scheme parameters, are indicated (the reader 
is referred to the original articles for details). 

6. Exact flow . 
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Approach and approximation 


V 


V 


w 


Wegner-Houghton eq. 113, 0(9°), numerically 





0.687 


0.595 


Wegner-Houghton eq. E3, 0(5°), polynom., M = 7 





0.657 


0.705 


Eq. for the Legendre effective action Eil 
0(9°), numerically 
0(9^), numerically 




0.054 


0.660 

0.618 


0.628 
0.897 


Polchinski eqE3 
0(9°), numerically 
0(9'^), numerically 




0.019-0.056 


0.649 

0.616 - 0.637 


0.66 

0.70 - 0.85 


World best estimates 


0.035(3) 


0.631 (2) 


0.80(4) 



Table 1: Results of calculations of the anomalous dimension i] and the critical exponents v and w for the Wilson-Fisher 

fixed point at d = 3 by various authors fm different versions of the ERG equation. The entries of the-ilast row (taken from 

the article by MorrisEJ'Cj) were obtained by averaging the world best estimatesE3. 



^^M{i)4- 



2M 



(34) 



Critical exponents characterize the flow very close to a particular fixed point. Within the approxima- 
tion, considered here, one can also study the flow globally. Thus, the Wegner-Houghton flow equation was 
studied in the leading order of the derivative expansion with the potential V{(j)^t) being approximated 
by the polynomial 

yM{(t), t) = -Ci(i)02 ^ _c2(t)04 + . 

(cf. (pl|))E3. We refer the reader to this article for details. 

7. Zaniolodchikov c-function . 

Another interesting application of the ERG equation is the calculation of the approximate Zamolod- 
chikov c-function which characterizes the geometry of the space of interactions of a given system. 

The Zamolodchikov c-function C{t) is a function which decreases monotonically along the flow (i.e. 
as the flow parameter t grows) and which is stationary at fixed points of the theory: 



dC{t) 
dt 



I fixed 



point 



0. 



Zamolodchikov proved the existepce of such function, which is unique up to a multiplicative factor, in 
two-dimensional unitary theories c3. There were some attempts to prove the existence of such function 
or construct it perturbatively in other cases cJ. Here we will show that for the Wegner-Houghton ERG 
equation the flow for the potential, approximated by the polynomial (p3), is gradient and permits a 
c-function description. The former means that the beta- functions of the couplings Ci{t), parametrizing 
the potential, are gradients of some function C(ci, C2, . . .): 



d , ^ dC 

-Q(t) = -.9., — , 



dt 



(35) 



where gij is a positive definite metpr; in the space of coupling constants. This implies that the set of 
renormalization flows is irreversible t£l. Then C is the c-function: 



dC 
'dt 



dC dci 
dcj dt 



<0. 



The potential 1^(0, t) in the polynomial approximation (|34| ) with M = 2 for the Wegoi 
equation was analyzed and the flow was shown to be gradient in this approximationtj 
system of flow equations consists of two equations: 



!r-Houghton ERG 
. In this case the 



Ci(t) 
C2(t) 



2ci 



6C2 

1 + Cl 



(4 - d)c2 



18ci 

(1+^ 



and can be cast into the form (Esl) with 



1 /I 



and 



^^■'" (TT^ I ^c 



± . ^ \/l " I ^ N*? ^ ! ^ \9 ^' ^9 



C(ci, C2) = --(1 + cif +3(1 + ci)3 - 3C2(1 + ci)2 + ^. (36) 



The c-function (jS^) is stationary at the fixed points. For example, for d = 3 it has the maximum at the 
Gaussian fixed point (ci = 0, C2 = 0) with C(0,0) = 1/6 and a saddle point at the Wilson-Fisher fixed 
point (c\ = -Ijl.cl = 2/49) with C{c\,cl) = 252/2401 w 0.105. 

5 Fermionic theory: first attempts. 

In this-secjtion we discuss briefly the results of ERG studies carried out for a two-dimensional fermionic 
model EaO. 

The Polchinski version of the ERG equation for pure fermionic theories can be derived in the same 
way as the one for the scalar theories (with another functional identity, of course). We start with a 
general action given by 

5 = y d''v^{v)Pl\v)^{-V) + 5™t(*, *, A), (37) 

with the regularized propagator equal to 

The analog of Eq. (ph is: 



dS_ 
'dt 



d'^q{2TTYK'{q 



dT^I(„2\ 



SS ^ SS 6 ^ 5S 



H{q) H{-q) Si){q) 6Tp{-q)_ 

, o f ,d fl-d + V ^ 2K'{q^\ fj, , SS ,, , SS \ 

''^'(*<'>'''^ji)-*<'""^5i))- <^«> 

Similar to the scalar case, rj is the anomalous dimension, ^p and ip are dimensionless fermionic fields, q is 
the dimensionless momentum and t is the renormalization flow parameter t = — ln(A/Ao). 

A concrete theory, studied within this approach, is the two-dimensional Euclidean chiral Gross- 
Neveau type model. This is a theory of fermions with N flavours described by 2-component spinors 
tp'^iq), '4'°'{q), a — 1,2, ... ,N, with SU{N) symmetry with respect to flavour indices of the left and right 
fields. The 7-matrices are given by Pauli matrices 71 = n, 72 = r2, and t^ plays the role of the chiral 
matrix 75. It can be shown that the most general action respecting such symmetries can be constructed 
out of the following operators: 

Siqi,q2) = r{qi)r{q2), P{qi,q2) = r{qi)l^r{q2), ^^(gi, 92) = ^"(gi)7^^"(92). (39) 

The ERG equation is solved for the truncated action which is represented by a finite sum S = S**-^'^^ + 
5'(4,o) _|_ 5'(4,2) ^ _ ^ 5'(n,m) ^^jj which contains all possible operators with up to n fermionic fields and 
up to m derivatives. Such expression can be viewed as a certain order of the derivative expansion with 
each term being further approximated by a polynomial in powers of fields. The general form of the first 
few S'("'™) respecting the symmetries of the theory can be shown to be 

S^^''^ = - j ^i^{-q)iqH<i). (40) 

^(4,2) - / '^'%'^f'^^ i'^^fKY. g») [gi(0(^(gl, g2)5(g3, g4) - P(gl, g2)P(g3, ^4)) 
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+ 52(^)^^(91,92)^^(93,94)], 

^(4,2) ^ J '^''^j^^^f g^ {2nf5{Y^ q,) [(mi (0(91 + 92)^ + m,{t){q, + q,){q, q,) 
+ TO3(i)(gi -92)^) {S{qi,q2)S{q3,q4,) - P(9i, 92)^(93, 94)) + ...]. 
The coefficients gi{t), mi{t), etc. are running coupling constants of various operators. Witfiin such 



approximation the ERG equation (38) becomes a system of equations for these couphng constants. S'^-^'^^ 
is of course the first term of the expansion of the kinetic part of ( p7| ) in powers of momenta (as before, 
we use normafization K{0) — 1). We assume that the normahzation of the kinetic term is fixed to be the 
canonical one, so no i-dependent coefficient appears here. 

The action with just S^^'^^ and S^^'^^ terms of ( [io| ) is the chiral Gross-Neveau modelE3. This model 
was studied in a number of papers and was shown to have non-trhqal fixed points which cannot be 
obtained neither by perturbative methods nor in the 1/A'^ expansion E3. One fixed point corresponds to 
the abelian Thirring model with gi = and 32 arbitrary. For the second one the coupling gi — 47r/(A^+ 1) 
and 92 is again arbitrary. 

It can be shown that the lowest order approximation which gives a fixed point solution with non-zero 
value for the anomalous dimension must contain at least 6 fermions and 3 derivatives. Execisely such 
action S = S'^^.i) + ^(^.s) + S^^'"^ + S'^'''^) + s-Cea) + 5(6,3) ^^^g ^^^^ ^j. calculations in Ref. N. It involves 
107 operators in total: 1 with 2 fermions and 1 derivative (the kinetic term), 1 in 5^^'^', 2 in S''^'^\ 11 in 
5^"*'^', 5 in 5^^'^^ and 87 in 5^^'^^. It turns out that the coupling constant of the operator with 2 fermions 
and 3 derivatives decouples from the rest of the system and does not play any role. With this action 
fixed point solutions were found and corresponding critical exponents were calculated. The results are as 
follows. 

First of all, consistent limits A^ ^ 00 of the system of the fixed point equations can be considered. 
In this limit the system simplifies dramatically. Two different large- TV regimes of the behaviour of the 
coupling constants were found and correspondingly one fixed point solution for each regime was obtained 
(let us call them fixed points I and II). 

Second, for finite A'^ > 1 among numerous solutions (most of which are spurious) two sequences of 
fixed points in the space of coupling constants were identified. Namely, as N grows these solutions match 
fixed point I and II. For each fixed point, both for finite TV and for N — 00, the corresponding anomalous 
dimension rj and the critical exponent Ai of the most relevant operator were calculated. The fixed point 
solutions are represented by their values of 77 and A in Fig. 1 and Fig. 2. 

For the sequence of solutions, which approach fixed point I, 77 ~ 1/TV, whereas Ai approaches the 
finite limit Ai — 1.231... (Fig. 1). The value of Ai for A^ = 00 is scheme independent for a wide range 
of scheme parameters. For this fixed point the constapts 172 and mi remain unfixed, a feature which 
resembles the exact result by Dashen and FrishmanEj. For finite N the value of Ai depends on the 
scheme, so the calculation was done for the regulating function K{z) = exp(— z^). 

The second sequence of solutions, the one which matches fixed point II, corresponds to the upper 
branches in Fig. 2. It exists only for N > 143. At the value N = 142.8 it joins another line of fixed point 
solutions. The values of 77 and Ai are scheme dependent and were fixed in some cases by the minimal 
sensitivity criterion. All coupling constants are fixed for these fixed points. 

The case iV = 1 was considered separately. In this case due to the Fierz relations the number of 
independent operators reduces. For example, all operators in 5^^'^' vanish because of the Grasmannian 
nature of the fermionic fields. In this case the approximation under consideration becomes the complete 
derivative expansion with terms up to 3 derivatives, i.e. there are no other operators with 1, 2 and 3 
derivatives apart from the ones described above. In this case the system of the ERG equations has just 
one fixed point solution. Values of all constants are fixed for this solution. 

Let us note that in this respect the case of finite iV > 1 is qualitatively different: in this case the 
approximation considered here is a truncation of the derivative expansion in the number of fields. In this 
case one observes a big number of fixed point solutions for each N. Most of them are, of course, spurious, 
similar to the polynomial approximation in the scalar case. Since successive approximations were not 
studied, it was not possible to app 
as it was done in the scalar theory 
matching type I or type II solutions when N —^ 00 were chosen 



the criterion of stability of solutions for various approximations, 
*^^' So, only those fixed points which form a sequence of solutions 
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6 Conclusions 

We would like to finish with a few conclusive remarks. 

1. The RG functional equations of type O) and the ERG equation (O) both reflect functional 
self-similarity of corresponding quantities and, thus, have very much in common at the fundamental 
level. The derivation of the RG equation essentially relies on the multiplicative character of finite Dyson 
transformations. They appear because of the freedom in fixing of finite arbitrariness which is left after 
the removal of ultraviolet divergencies. This means, first, that some underlying perturbative expansion in 
a coupling constant (or analyticity in the coupling constant) is assumed and, second, that by construction 
the RG approach is formulated for renormalizable theories with the upper cutoff Aq being sent to infinity 
and, thus, all irrelevant operators being removed. This is not the case for ERG approach where the cutoff 
is kept finite. Thus, the built-in functional self-similarity of functional RG equations of type (]^ imposes 
strong restrictions and reduce the number of arguments of g{xjy.g), but these equations do not contain 
enough of the dynamical information to determine the functional. The only regular way to proceed is to 
use renormalization group functions (Gell-Mann-Low /3-function, etc.) calculated within the perturbation 
theory where contributions of irrelevant operators enter through loop corrections. Contrary to this, the 
ERG equation ( |l3| ) basically contains information about renormalization group evolution of the whole 
Wilson effective action, i.e. about all operators and, as such, does not need any additional inputs. This 
allows to search for non-perturbative (i.e. non-analytical in the coupling constant) solutions Ij. 

2. The ERG method was demonstrated to be a powerful approach for non-perturbative studies of the 
continuum limit in quantum field theory for scalar and fermionic theories. 

3. The proper derivative expansion is an effective and reliable technique which allows to search for fixed 
point solutions and calculate critical exponents. When combined with further polynomial approximation, 
some qualitative features can be reproduced. However, numerous spurious solutions appear and a special 
procedure should be applied to identify the true fixed points. 

4. For ERG equations with an arbitrary cutoff function K{z) results beyond the leading 0(9°) ap- 
proximation depend on K{z). This regulator dependence is similar to the scheme dependence in usual 
perturbation theory in the RG approach. This is also related to an important issue of reparametrization 
invariance. We_did not consider it here, and the readers are referred to other articles where this problem 
is discussed Ejo. 

5. In spite of some progress in the ERG approach in gauge theories cHI, further developments should be 
made before we have a regular tool for obtaining non-perturbative quantitative results for this class of 
theories. 
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Figure 1: The values of A^r; (solid line) and Ai (dashed line) for the sequence of fixed points approaching type I fixed point 

of the N = oo case. 
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Figure 2: The values of rj (solid line) and Ai (dashed line) for the second sequence of fixed points, 
approach the corresponding values of type II fixed point of the N = oo case. 



The upper branches 
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